Orbital order in vanadium spinels 
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Motivated by recent theoretical and experimental controversy, we present a theoretical study to 
clarify the orbital symmetry of the ground state of vanadium spinel oxides AV2O4 (A=Zn, Mg, 
Cd). The study is based on an effective Hamiltonian with spin-orbital superexchange interaction 
and a local spin-orbit coupling term. We construct a classical phase-diagram and prove the complex 
orbital nature of the ground state. Remarkably, with our new analysis we predict correctly also 
the coherent tetragonal flattening of oxygen octahedra. Finally, through analytical considerations 
as well as numerical ab-initio simulations, we propose how to detect the predicted complex orbital 
ordering through vanadium K edge resonant x-ray scattering. 

PACS numbers: 75.10.Jm, 75.30.Et 



Vanadium and titanium spinels, AB2O4 (B=Ti 3+ , or 
V 3+ ), belong to a class of frustrated antiferromagnets 
where magnetic B-ions are characterized by orbital de- 
generacy due to partial occupancy of t2 S -orbitals (n t2 =1 
for titanates and n t2g =2 for vanadates). Recently these 
spinels were thoroughly studied from both experimental 
000 and theoretical 0, points of view. While 
the ground state of Ti-based spinels can be explained in 
terms of orbitally-driven superexchange interactions on 
the frustrated pyrochlore lattice 0, the situation seems 
not so fluid for vanadium spinels, as two conflicting the- 
oretical works appeared to explain their structural and 
magnetic properties . 

In AV2O4, magnetically active V 3+ -ions form a py- 
rochlore lattice and are characterized by two 3d electrons 
in t2 ff -orbitals, while A is a divalent ion like Cd 2+ , or 
Zn 2+ , or Mg 2+ . All compounds show qualitatively sim- 
ilar structural and magnetic behavior with a structural 
transition at a higher temperature T$ and an antifer- 
romagnetic (AFM) transition at a slightly lower tem- 
perature Tjv 0. These findings have been interpreted 
by Tsunetsugu and Motome [J] as an interplay of dda 
superexchange (SE) interaction and geometrical frustra- 
tion: they showed that ordering of orbitals can partially 
remove magnetic frustration and explain the experimen- 
tally observed magnetic structure which is composed of 
AFM chains running in [110] and [110] directions. The 
ground state orbital ordering suggested in Ref. con- 
sists of stacked ab planes with alternating d xz and d yz 
vanadium hole orbitals (hereafter referred to as ROO). 

On the other side Tchernyshyov [3 pointed out that 
the ground state symmetry I4i/a of ROO solution seems 
at odds with x-ray and neutron diffraction data, indi- 
cating a I4i/amd space symmetry. Thus, he proposed a 
purely ionic model where spin-orbit (SO) coupling plays 
the major role and the V hole occupies (predetermining 



the sign of Jahn- Teller (JT) distortion) a complex linear 
combination of x z and yz orbitals: (d xz Hd yz )/V2. (We 
shall refer to this orbital order as COO). 

Actually, the correct space group of the system is still 
elusive. The tetragonal I4i/amd space group was found 
in Ref. [l|, while the authors of the neutron scattering 
experiment supported ROO (thus, I4i/a space group ) 
as more compatible with their findings. They found a 
clear signature of the one-dimensional character of spin- 
fluctuations in the temperature range T/v < T < Tg, 
indicating a weak coupling between AFM chains. Based 
on this observation, and claiming that interchain cou- 
pling is much stronger for COO compared to ROO, they 
interpreted the experimental results in favor of the latter. 

The aim of the present paper is to shed some light 
on this issue. To this end, we first link the two appar- 
ently unrelated pictures emerging from Refs. 00 and set 
them in an unique framework, by considering SE interac- 
tion and SO coupling on an equal footing and construct- 
ing a classical ground state phase diagram. We demon- 
strate that the SO coupling favors the states with un- 
quenched orbital momentum for any value of the coupling 
strength. We show that the only ground-state which is 
compatible with the experimentally observed magnetic 
structure is the COO-phase. This state is also charac- 
terized by a coherent tetragonal flattening of the oxygen 
octahedra along the same axis, as in real materials. Re- 
markably, this result is not pre-determined by any ad hoc 
choice of JT interaction as was done in Refs. [4j,la], but 
simply follows from our treatment of SE and SO interac- 
tions. We also evaluate the interchain spin couplings for 
both ROO and COO and find that they are comparable. 
Thus the experimental results of Ref. [2| cannot discrim- 
inate between the two types of orbital ordering. Finally, 
we suggest a possible key experiment to single out the 
proposed complex orbital states by means of magnetic 
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resonant x-ray scattering (RXS). To this aim, we per- 
formed a formal analysis and an ab-initio numerical sim- 
ulation by means of the relativistic multiple scattering 
code in the FDMNES package 0,0. 

The effective Hamiltonian for AV2O4. — The model 
Hamiltonian consists of two parts H = Hse + Hso- The 
first term represents the superexchange interaction be- 
tween vanadium ions and the second term stands for 
the local spin-orbit coupling. First we discuss Hse- 
The insulating phase of AV2O4 is of Mott-Hubbard type 
and the SE interaction can be described by the Kugel- 
Khomskii model We borrow the SE Hamiltonian 

for V 3+ ions, with three-fold degenerate t^g orbitals and 
spin S=l, from Ref. ^lj, to which we refer for tech- 
nical aspects, and specialize it to our case. Hund's 
exchange, Jh — 0.68 eV, and Coulomb on-site repul- 
sion (same orbital), U\ ~ 6.0 eV, are taken from spec- 
troscopy data |l_2j . For future convenience, we intro- 
duce the small parameter 77 = Jh/U\ ~ 0.11. Tight- 
binding fit to linear-augmented plane-wave calculations 
[l3| shows that dominant overlap integrals are of cr-type 
with t = i/U dda = -0.24 eV For this reason, we 

consider only nearest-neighbor (NN) hopping integrals of 
dda-kind on the pyrochlore lattice, as in Refs. 4, 6]. The 
dda overlap in a/3 plane (af3=xz, yz, xy) connects only 
the corresponding orbitals of the same a/3 type. Thus, 
orbital operators are Ising-like and their contribution can 
be expressed simply in terms of projectors Pi ja /3 onto the 
occupied orbital state a/3 at site i, which are equal to 1 
when a/3 orbital is filled, or to when it is empty. We 
finally get: H S e = § £V Hy, where 

Ha = - [J Si ■ Sj + Ji] On -J 2 [i-Si- Sj] da (1) 

is the energy of bond ij. The orbital contributions 
Oij and Oij along the bond ij in a/3-plane are given 
by: 0^ = Pi, a p(l - Pj, a p) + Pj,ap(l - Pi,ap) and 
0^ = Pi, a pPj, a /3- We defined: J = r?J/[l - 3rj\, 
Ji = J[l - n]/[l - 3r7], J 2 = J[l + n]/[l + 2??], and the 
inequality Jo < J 2 < Ji holds. The overall energy scale 
is J = tl d jU x ~ 9.6 meV. 

The remaining local part of the effective Hamiltonian 
describes the relativistic spin-orbit coupling at each vana- 
dium ion. We adopt the usual representation of effective 
V = 1 orbital angular momentum for ^-electrons |l 5| . 
with L' z \xy) = 0, and L' z \xz±iyz) /y/2 — ±\xz±iyz) / ' \[2. 
The true angular momentum L is related to the effective 
one by L ~ —L'. With this notation Hso can be written 
as: H so = -A £ 4 1\ ■ S h with A > 0. 

Real orbital configuration. — We first consider the 
states when only real orbitals are occupied. The orbital 
angular momentum is quenched and the local spin-orbit 
coupling is inactive. Such orbital patterns are eigenstates 
of the orbital part of Hse- This allows to classify all pos- 
sible interacting bonds ij in a/3 planes in two categories: 




FIG. 1: Orbital, spin and bond arrangements on a tetrahe- 
dron, (a), (b), (c), (d), and (e) respectively correspond to 
ROO, ROO-I, COO-I, COO-II, and COO phases discussed 
in the text. Hole orbitals are shown. The solid (dashed) 
lines stand for AFM (FM) bonds. The thicker is the line the 
stronger is the coupling on the bond. 

i) weakly ferromagnetic (FM) 61-bond, when only one of 
the two sites has an a/3 electron {Oij contribution); and 

ii) strongly AFM & 2 -bond with both i and j ions occupied 
by a/3 electrons {Oij contribution). There are also non- 
interacting bonds when nor i neither j have a a/3 orbital 
occupied. As we are interested in magnetic configurations 
of the whole pyrochlore lattice, in the following we limit 
our considerations to classical spins. The exchange cou- 
pling on 61-bonds is given by H^ = —JoSi ■ Sj — Ji, with 
maximum FM energy gain: — Jq + Ji . For r) = the 
bond is nonmagnetic and = J. The exchange cou- 
pling of 62-bonds is given by Hb 2 = — J 2 (l — Si ■ Sj), and 
the maximum (AFM) energy gain is equal to Eb 2 = 2J 2 . 
For future convenience we point out that E{, 2 = 2£'b 1 for 
i] = 0, while Eb 2 < 2E) >1 for any finite r\. 

In order to identify the lowest energy state in the lat- 
tice, we introduce the operators iV^ = Y^,{ij) Oij and 
^b 2 ~ S(-y) Oij that count the number of interacting 61 
and 62 bonds on a tetrahedron. Using the identity for 
projectors Pi >X y + Pi,xz + Pi,yz — 2 (as there are two elec- 
trons occupying different orbitals at each site), one ob- 
tains + 2N b T — 8. Rewriting this constraint in terms 
of energies of interacting b\ and 6 2 bonds, one finds that 
the energy gain per tetrahedron for 77 = is always equal 
to JN^ + 2JN b T > = 8 J . Therefore, for r) = 0, all real 
orbital patterns are degenerate (lfij |. However, as with 
increasing 77 the energy gain on & 2 -bond decreases, while 
that on b\ increases, the degeneracy is lifted. Thus, the 
lowest energy configuration is the one with the maximum 
number of b\ bonds per tetrahedron, which is = 4 
(recall that + <6-the total number of bonds on 
a tetrahedron). 

There are two topologically different tetrahedral con- 
figurations with = 4 [see Fig. [Ja), Hfb)]. In the 
case of ROO (Fig. da)) there are two strongly AFM 
& 2 -bonds located on opposite edges and coupled through 
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FIG. 2: Classical phase diagram of Hse + Hso in (A, rf) plane. 
The shaded rectangle denotes the domain of possible values 
of parameters. 

weakly FM (J' = — Jo), frustrated &i-bond interactions. 
When all tetrahedra are of ROO-type, the average energy 
per site is £roo = — 2Ji — 2J2. The magnetic structure 
can be viewed as a collection of AFM chains running in 
[110] and [110] directions, coupled through frustrated in- 
terchain interaction J'. At the classical level the chains 
are decoupled from each other. 

In the configuration ROO-I the two strong AFM bonds 
share a corner and not all 6i-bonds are frustrated (Fig. 
H^b)). The spins of each tetrahedron form a fully collinear 
up-up-up-down (uuud) state with a finite magnetic mo- 
ment equal to 1/2 of saturation value. The energy of this 
state -Eroo-i = — Jo — 2Ji — Ui is the lowest among 
those with quenched angular momentum. Thus, ROO-I 
phase is the ground-state solution when A = 0. How- 
ever, it is at odds with experiments in both magnetic 
structure and JT distortion. In fact, the latter should 
have the same local axis at all sites, while this is not the 
case for ROO-I phase. Notice again that in Ref. 4] ROO 
ground state was found because JT coupling with only 
tetragonal distortion in z-direction was pre-introduced. 
However, the inclusion of a local JT coupling with cubic 
symmetry will not affect the above results, as each state 
from the manifold would gain the same amount of JT 
energy by an elongation of each V06 octahedron along 
the local tetragonal axis [Io| . 

Complex orbital configuration. — Now consider the or- 
bital configuration shown in Fig. 1(c) (COO-I), obtained 
from ROO-I by replacing a yz hole with [xz + iyz}/^/2 at 
site 1. This site now carries angular momentum L' z = 1 
and gains the SO energy A. The spin couplings on bonds 
1-2 and 1-4 arc obtained by replacing orbital bond op- 
erators Oi2(4) and Oi2(4) by their expectation values 
(012(4)) = (Oi2( 4 )) = 1/2. This leads to the spin SE 
Hamiltonian H 12 (4) = [- Ji - Ji + (J2 - Jq)Si • S 2 (4)]/2, 
with AFM coupling [J2 — Jo]/2 > 0. Its ground-state en- 
ergy is i?coo-i = — Jo — 2 J\ — 2J 2 — A/4, which is lower 



than the lowest energy for real orbital patterns (ROO- 
I discussed above). Therefore, for any finite value of A 
the states with quenched angular momentum are unsta- 
ble with respect to those with unquenched L' z . We can 
further construct the states with complex orbitals clas- 
sifying them by the number of sites with non-zero L' z 
per tetrahedron. The state with two complex sites (state 
COO-II, shown in Fig. 2(d)) has an average energy per 
site given by -Ecoo-11 = — 3J2 — J\ — A/2. It is lower 
than the energy of COO-I for A > 4[Ji - J 2 + J ] ~ 877. 
The magnetic order for such an OO can be regraded as a 
quasi-two-dimensional AFM pattern, formed by strong- 
bonds (see the path 1-2-4-3-1 in Fig. 1(d)), with weak 
FM interaction along diagonals (bonds 1-4 and 2-3). The 
configuration with three complex sites has never the low- 
est energy. When all four sites carry non-zero L' z (state 
COO, shown in Fig. 1(e)), the ground state energy is 
Ecoo — — 1/2[5J2 + 2Ji] — A. The tetrahedron is char- 
acterized by two strong AFM 62-bonds (2-4 and 1-3) 
and the lattice magnetic structure is given by ID AFM 
chains running in [110] and [110] directions. The in- 
terchain coupling is frustrated and weakly AFM, with 
J' = 1/4 [J 2 — 2 Jo], contrary to the ROO-phase whose 
J' = —Jo is FM. Its strength is practically the same 
for both COO and ROO. For example, for r\ = 0.12, 
I J'| = 0.14J for ROO and | J'\ = 0.16J for COO. There- 
fore, in both cases the magnetic network is quasi-onc- 
dimensional and the experimental results of Ref. can 
not rule out one or the other, as we have already out- 
lined. Notice that the frustration of COO-phase is much 
weaker than ROO phases, as the spin direction is fixed 
along the z-axis by SO coupling and the degeneracy of 
the ground state is discrete as discussed in Ref. [j| . 

In Fig. [21 the ground-state phase diagram is shown in 
terms of the two dimensionless parameters A = A/ J and 
rj. The shaded rectangle represents the domain of real- 
istic values for parameters: 77 = 0.8 — 0.12 0,^3; an d 
A = 13 - 20 mcV (A = 1.35 - 2.08) (3,113 Interestingly, 
if the system were close to the phase boundary, the ap- 
plication of an external magnetic field might induce the 
transition from COO, with zero magnetization, to COO- 
I, with a finite magnetic moment. The value of the field 
depends on how far from the phase boundary actual pa- 
rameters locate the system. For the parameters within 
the shaded rectangle in Fig. [5] the maximal value of the 
critical field is estimated to be H c ~ 86 T. The magnetic 
field experiment can be an interesting test of the phase 
diagram presented here. 

To conclude the above consideration, we note that 
among ground states appearing in the phase diagram, 
only COO is compatible with the experimentally ob- 
served AFM structure. This type of orbital ordering 
also implies the coherent flattening of V06 octahedra 
along a unique axis leading to the experimentally ob- 
served tetragonal distortion. In what follows, we discuss 
a key experiment that we believe can give unambiguous 
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way to detect the suggested ground state orbital order. 

Experimental detection of COO-phase. — An effective 
way to unravel the question of the orbital symmetry, both 
below TV and between Tjv and Tg, is by means of V K 
edge RXS, where one is sensitive to the on-site value of 
L z , as the local transition amplitudes are added with a 
phase factor that compensates the vanishing effect due 
to the global symmetry 0. The elastic scattering am- 
plitude is A(Q, w) = £j-e l< 2' Pj : fj{u), where Q is the mo- 
mentum transfer in the scattering process, hoj the pho- 
ton energy, pj the atomic positions in the unit cell, and 
the sum is over the 8 atoms in the tetragonal cell. The 
atomic scattering factors (ASF), fj{u>), is a second-order 
process in the electron-radiation interaction, give n by (in 
the dipole approximation and atomic units) : 
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FIG. 3: Magnetic RXS at (lll) CT7r reflection for CdV 2 4 . 
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where |\&o ) (|*n')) is the ground (excited) state, with 
the origin taken on the j-th scattering atom, and Eq {E n ) 
is its energy; r„ takes into account the finite lifetime of 
the excited states; £i( s ) is the incident (scattered) polar- 
izations, and r is the coordinate of the electron in the 
reference frame of the resonant ion. In COO-phase the 
space-group is I4i/amd, and we can write the structure 
factor in terms of the ASF of one single ion, through the 
symmetry operators C 2z and Cf z (two- and four-fold ro- 
tation around c-axis) and T (time-reversal): S{Q) = (l+ 

(_)*+*+«£■) [i + {-f+ l c 2z + i h - k+l c+ z + i h+k - l c£]h. 

For the I4i/a space group, induced by the reduced 
symmetry of the ROO-phase, we get: S{Q) = (l + 
{-) h+k + l f)[l+{-) h + l C 2z +i- h+k+l C++i- h - k - l C^]f 1 . 
Therefore, for h+k+l=odd, both signals are proportional 
to (1 — T), As K edge magnetic signals are a signature 
of non-zero orbital angular momentum, we can use V K 
edge RXS to discriminate between ROO phase, where no 
signal is expected, and COO phase, where the signal is 
proportional to L z [9j- 

We have then performed a numerical simulation to es- 
timate the order of magnitude of this signal for complex 
orbital occupancy with the relativistic multiple scattering 
code in the FDMNES package Results are shown in 
Fig. [3]for the (111) reflection in CdV204 (structural data 
were taken from Ref. Q). The incoming polarization is 
directed along one of a, b axes, where the intensity is the 
highest. The signal is of the order of 10 _3 rQ at dipolar- 
energies (ro is the classical electron radius) , which is well 
within the present capabilities of third generation x-ray 
sources: this is a key experiment to choose the ground 
state in vanadate spinels. Probably a detectable signal 
could be obtained also through non-resonant magnetic x- 
ray scattering, thus allowing a quantitative measurement 
of L z . 
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